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Abstract 



We study a systematic improvement of perturbation theory for gauge fields 
on the lattice; the improvement entails resumming, to all orders in the cou- 
pling constant, a dominant subclass of tadpole diagrams. 

This method, originally proposed for the Wilson gluon action [1], is ex- 
tended here to encompass all possible gluon actions made of closed Wilson 
loops; any fermion action can be employed as well. The effect of resummation 
is to replace various parameters in the action (coupling constant, Symanzik 
coefficients, clover coefficient) by "dressed" values; the latter are solutions to 
certain coupled integral equations, which are easy to solve numerically. 

Some positive features of this method are: a) It is gauge invariant, b) it 
can be systematically applied to improve (to all orders) results obtained at 
any given order in perturbation theory, c) it does indeed absorb in the dressed 
parameters the bulk of tadpole contributions. 

Two different applications are presented: The additive renormalization of 
fermion masses, and the multiplicative renormalization Zy {Za) of the vector 
(axial) current. In many cases where non-perturbative estimates of renor- 
malization functions are also available for comparison, the agreement with 
improved perturbative results is significantly better as compared to results 
from bare perturbation theory. 

Keywords: Lattice QCD, Perturbation theory. Improved actions. Tad- 
pole improvement. 
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I. INTRODUCTION 



Since the earliest studies of quantum field theories on a lattice, it was recognized that 
quantities measured through numerical simulation are characterized by significant renormal- 
ization effects, which must be properly taken into account before meaningful comparisons 
to corresponding physical obscrvablcs can be made. 

As has been rigorously demonstrated [2] , the renormalization procedure can be formally 
carried out in a systematic way to any given order in perturbation theory. However, calcu- 
lations are notoriously difficult, as compared to continuum regularization schemes; further- 
more, the convergence rate of the resulting asymptotic series is often unsatisfactory. 

A number of approaches have been pursued in order to improve the behaviour of per- 
turbation theory, among them Rcfs. [3,4]. These approaches share in common the aim to 
reorganize perturbative series in terms of an expansion coefficient which would be more 
suitable than the bare coupling constant go; the definition of such a "renormalized" cou- 
pling constant is not unique, but can depend on the observables under study and on an 
energy scale. It is expected that such a definition will reabsorb a large part of the tadpole 
contibutions which are known to dominate lattice perturbation theory. 

Some years ago, a method was proposed to sum up a whole subclass of tadpole diagrams, 
dubbed "cactus" diagrams, to all orders in perturbation theory [1,5]; this procedure has a 
number of desirable features: It is gauge invariant, it can be systematically applied to 
improve (to all orders) results obtained at any given order in perturbation theory, and it 
does indeed absorb the bulk of tadpole contributions into an intricate redefinition of the 
coupling constant; in cases where non-perturbative estimates of renormalization coefficients 
are also available for comparison, the agreement with cactus improved perturbative results 
is significantly better as compared to results from bare perturbation theory. 

In the present work we extend the improved perturbation theory method of Refs. [1,5], to 
encompass the large class of actions which are used nowadays in simulations of QCD. This 
class includes Symanzik improved gluon actions with any arbitrary combination of closed 
Wilson loops, combined with any fermionic action. In Section II we present our calculation, 
deriving expressions for a dressed gluon propagator, as well as for dressed gluon and fermion 
vertices, as a result of the summation of cactus diagrams to all orders. We show how 
these dressed constituents are employed to improve 1-loop and 2-loop Feynman diagrams 
coming from bare perturbation theory. In Section III we apply our improved rcnormalizaton 
procedure to a number of test cases involving Symanzik gluons and Wilson/clover/overlap 
fermions. Improvement of QED is relegated to an Appendix. 

Clearly, all resummation procedures, whether in the continuum or on the lattice, bear 
a caveat: A one-sided resummation could ruin desirable partial cancellations which might 
exist among those diagrams which are resummed and others which are not; what is worse, 
the end result might depend on the gauge. As we shall see, no partial cancellations will be 
ruined in our procedure, due to the distinct A^-dependence of the resummed diagrams {N is 
the number of colors); furthermore, our results will be gauge independent. 
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II. THE METHOD 



In this Section, following the outline of Ref. [1], we start illustrating our method by 
showing how the gluon propagator is dressed by the inclusion of cactus diagrams. We will 
then dress gluon and fermion vertices as well. Finally, we will explain how this procedure 
is applied to Feynman diagrams at a given order in bare perturbation theory, concentrating 
on the 1- and 2-loop case. 



A. Dressing the propagator 



We consider, for the sake of definiteness, the Symanzik improved gluon action involving 
Wilson loops with up to 6 links; it will be evident from what follows that the method 
is applicable to any gluon action made of Wilson loops. In standard notation (see, e.g., 
Ref. [6]), the action reads: 



. 2 

'G — ^ 

90 



Co ReTr (1 - t/piaquette) + Ci ^ Re Tr (1 - t/rectangle) 

plaquette rectangle 



+ C2 ^ Re Tr (1 - t/chair) + C3 Re Tr (1 - t/parallelogram) 

chair parallelogram 



The coefficients q can in principle be chosen arbitrarily, subject to a normalization condition 
which ensures the correct classical continuum limit of the action: 



Co + 8ci + 16c2 + 8C3 = 1 (2) 

Some popular choices of values for Cj used in numerical simulations will be considered in 
the apphcations of Section III. The quantities C/j {i — 0,1,2,3, respectively: plaquette, 
rectangle, chair, parallelogram) in Eq. (1) are products of link variables Ux,n around the 
perimeter of the closed loop. 

Applying the usual parameterization of links in terms of the continuum gauge fields 

C/,,^ = exp(i5oaA^(x + a/i/2)), = T", Tr (T" T") = | (3) 

(a : lattice spacing, set to one from now on; fi : unit vector in direction fi; T'^ : generator of 
SU{N) algebra) and the Baker- Campbell-Hausdorff (BCH) formula, Ui takes the form: 

U, = exp (i ^0 fI'^ + i gl + i + 0{g',)) (4) 

where Fj;^^ is simply the sum of the gauge fields on the links of loop i (e.g., for the plaquette: 
F^^^ = A^{x+il/2) + A^{x+fi+u/2) - A^{x+u+il/2) - A^{x+u/2)), while F^^^^ (j > 1) are 
j-th degree polynomials in the gauge fields, constructed from nested commutators. 
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We may now define the cactus diagrams 
which dress the gluon propagator: These are 
tadpole diagrams which become disconnected 
if any one of their vertices is removed (see Fig. 
1); further, each vertex is constructed solely 
from the f/^'* parts of the action. 

A diagrammatic equation for the dressed 
gluon propagator (thick line) in terms of the 
bare propagator (thin line) and 1-particle ir- 
reducible (IPI) vertices (solid circle) reads: 

The IPI vertex obeys the following recursive equation: 




Figure 1: A cactus 
- + ■ ■■ 
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: (6) 

In order to put this equation into a mathematical form and solve it, let us first write down 
the bare inverse propagator resulting from the action (1), and from the gauge fixing 
term: 



S,f = Y^,E Tr {A;A^{x+fi/2) A;A,(x+z>/2)), A; 

> X,IJ," 



(x) = (f){x—fi) — 0(x) (7) 



The quadratic part of the total gluon action thus becomes (in the notation of Ref. [6]): 

1 r d^k 



5, 



(0) 



2 J-^ i2-ny 



T.^l(k)D-l{k)Al{-k) 



(8) 



where: Dj,{k) = ^ {kjS^^ - k^kpdp^^ d^p + kpk^ 

and: dp, = (1 - 5p,) [Cq - Ci P - {Pp + H 



kn 



ku = 2 sin — 
M 2 



The coefficients Co, Ci, C2 are related to the Symanzik coefficients q by 
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Co = Co + 8ci + 16C2 + 8C3 , Ci = C2 + C3 , C2 = Ci - C2 - C3 



(9) 



The inverse propagator can thus be put in the form: 

D;^{k) ^ Co + ci + C2 + C3 ^^(A;) + ^ k,K (10) 

The matrices G^*^ (k) are symmetric and transverse, i.e. they satisfy: 

EG»(fc)i = (11) 

Each of them originates from a corresponding term: Tr(F/^'* f/^'*) of the gluon action. Con- 
sequently, each of the diagrams on the r.h.s. of Eq. (6), being the result of a contraction with 
only two powers of f/^^ left uncontracted, will necessarily be equal to a linear combination 
of G^'^\k); this imphes that the IPI vertex G^^^{k) (the l.h.s. of Eq. (6)) can be written as: 

G^^\k) = ao G^^^ (k) + ai G'(^) (fc) + as G^^^ (k) + as G^^^ (k) (12) 

Each of the quantities will in general depend on N, go, cq, ci, C2, C3, but not on the 
momentum. We must now turn Eq. (6) into a set of 4 recursive equations for . 
Eq. (5) leads to the following expression for the dressed propagator D'^^{k): 

D''' = D + D G^Pi D + D G^Pi D G'''' D + --- = D (I3) 

=^ (L»dr-)-l ^ (j^ _ ^IPI ^-1 ^ ^-1 _ ^IPI 

= Co G^") + ci G'(^) + C2 G^') + C3 G(^) + ^ fc^i , Q = c, - a, (14) 

We observe that dressing affects entirely the transverse part of the inverse propagator, 
replacing the bare coefficients Cj with improved ones Cj , and leaves the longitudinal part 
intact. The same property carries over directly to the propagator itself; the consequence of 
this will be that our method leads to the same results in all covariant gauges. 

In terms of the dressed propagator, Eq. (6) can be drawn as: 

= SL + ^ + ^ + 

(15) 

Let us now evaluate a typical diagram on the r.h.s. of Eq. (15); it will be the sum of 4 
terms, one term for each of the Wilson loops Ui (i = 0, 1, 2, 3) in the action, from which its 
n-point vertex may have originated. There are {n — 2)/2 1-loop integrals in the diagram; 
each of them corresponds to the contraction of two powers of Pj^^^ via a dressed propagator, 
and will contribute one power of /?j(co, Ci, £2, C3), where: 
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/3: 



/32 = £ ^ (C(8-?1) D%{q)/2 - qM^-ql) D%{q)/2) 

Ps - £ ^ (3 ^Jp(?)/2 - 3 g, J(g)/2) (16) 

{fiji/jp assume distinct values; no summation implied). Once again, we note that Pi are 
gauge independent, since the longitudinal part cancels in the loop contraction. 

For the contraction of the SU (N) generators we first evaluate F{n; N), which is the sum 
over all complete pairwise contractions of Tr{T'"iT'"2 . . .T""}: 

F{2-N)=6a,a,Tr{T''^T''''} 

F(4; AT) = {da,a, 5a,a, + Kas Ka, + Ka, Kas) Tr {T^' T'^'^ T'^' T'^'} 



PeSn 



{F{2n + 1; A^) = 0; S'n is the permutation group of n objects). The generating function 
G{z; N) for this quantity: 

00 n Jn 

Giz; N)^Y: rT^(^; => F{n- N) = —G{z; N)\^^, (18) 

n=0 

has been computed explicitly in Ref. [1], using Gaussian integration over the space of trace- 
less Hermitian matrices [7], with the result: 

G{z; N) = e^'(^-i)/(^^) L],_,{-z'' /2) (19) 

{Lp{x) : Laguerre polynomials). Since 2 out of n generators remain uncontracted in our 
case, color contraction does not lead to F(n; N), but rather to: 

nF{n;N) 

2(iV2-l) ^^"^ 

Thus, upon contraction, an n-leg diagram in Eq. (15), with its vertex coming from the term 
Ui of the Langrangian (i = 0, 1, 2, 3), will merely result in the following multiple of G^'^^ : 

Q (iffo)" nF(n;N) (n-2)/2 ^(.) . 

We are finally in a position to set Eq. (15) in a mathematical form: 

G^'^ = 

^ ^ Q (^^?o)" nF(n; AT) 
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Unknown in Eq. (22) arc the coefficients Ui ; they appear on the Lh.s., as well as inside the 
integrals j3i of the r.h.s, by virtue of Eqs. (16, 14). We recall that G*-*^ arc functions of the 
external momentum k ; if these are independent^, then Eq. (22) amounts to 4 equations for 
the 4 coefficients . 

The generalization of our procedure for improved gluon actions with arbitrary numbers 
and types of Wilson loops is now evident. It is crucial to check at this stage that all 
combinatorial weights are correctly incorporated in Eq. (22); this is indeed the case. 

Splitting Eq. (22) into 4 separate equations, and making use of Eq. (18), we can recast 
the infinite summations in closed form: 

a^^ ~ 1 (^^o)" nF(n:iV) (,_2)/2 
a n=ks,...9i n\ 2(iV2-l) 

2 



(23) 



^ ^-^{N'-l) = ^G'{z;N) 

Ci Z 2=(iflo/3/ ) 

^ e-^-o (N-mm L],_M A/2) + 2 Ll_,{gl A/2)) (24) 

In solving Eqs. (24), each choice of values for (q , go , N) leads to a set of values for 
Ci = Ci — ai . The latter are no longer normalized in the sense of Eq. (2); one may equivalently 
choose, however, to express the results of our procedure in terms of a normalized set of 
improved coefficients, Ci/Co and an improved coupling constant = qI/Cq, where: Cq — 
Co + 8ci + 16c2 + 863 . In fact, it is convenient to treat bare and improved coefficients on an 
equal footing, by defining rescaled quantities as follows^: 

7i = ^, 7i=^, A(C0,C1,C2,C3) = ^0 A(C0,C1,C2,C3) = A(70,7l,72,73) (25) 

The rescaled quantities % must now satisfy the coupled equations: 

li - ^7.e-^'(^-^)/(^^) (^Lj,_,(A/2) + 2L^_,(A/2)) (26) 

For the gauge groups SU{2) and SU{2>), the Laguerre polynomials have a simple form, 
making Eqs. (26) more explicit: 



^Actually, G^'^^ is not independent of the rest, so that we have 3 equations for 3 coefficients; this 
causes no complication. In any case, typically C2 = in simulations. 

^The dressed propagators in A will now contain a rescaled gauge parameter (1—^) {1—^), 
which is irrelevant since the longitudinal part does not contribute. 
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(iV = 2) : 7. = 7. e-^^/' " 1) ' (^^ = 3) : 7. = 7. e'^^/' ' J + ^) (27) 

Given the highly nonhnear nature of Eqs. (26), it is not a priori clear that a solution for % 
always exists^; it turns out that this is always the case, for all physically interesting values 
of Q, and for all values of go ranging from (^o = up to a certain limit value, well inside the 
strong coupling region. 

Fortunately, numerical solutions of Eqs. (26, 27) can be found very easily. We use a fixed 
point procedure, applicable to equations of the type x = f{x) : 

7. = /.(7.) => 7. = Jhn^7t\ where: #=7., it^'^ = Wt^) (28) 

In order for the procedure to converge (attractive fixed point), it must be that: \dfi/d'yi\ < 1 
in a neighborhood of 7i. This has been verified in a number of extreme cases. 



B. Numerical values of improved coefficients 



Here we present the values of the dressed coefficients for several gluon actions of interest. 
I 1 1 1 1 1 1 1 1 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

Fig. 2: Improved coefficient cq for N=2 and A^=3 (plaquette action) 



• Let us start with the plaquette action (co=l, ci=C2=C3=0) [1]: In this case, Eqs. (26) 
reduce to only one equation, for 70, while 7i=7i=0 (z=l,2,3). This equation is further 
simplified greatly since the integral /5o can now be evaluated in closed form, Pq = 1/(270); 
for iV = 3 we obtain (cf. Eq. (27)): 

, -,oV(12£o) + (29) 



^The converse, of course, is trivial: Finding the bare values 7^ which lead to a given set of dressed 
values 7j is immediate, since the integrals Pi only depend on 7j, not 7^ . 
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In Fig. 2 we plot cq (in the notation of [1], cq = 1— w((?o)) as a function of (^g, for = 2 
and = 3. The range of Qq values, for which solutions exist, extends from gQ = (where 
Co = 1) up to 16y/e/3 ~ 3.23 (A^ = 2) and 1.558 (A^ = 3); this covers the whole region of 
physical interest. 

1.8 
1.6 
1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 

0.0 0.5 1.0 1.5 2.0 2.5 

Fig. 3: Improved coefficients cq and ci (tree-level Symanzik improved action) 

• The tree-level Symanzik improved action [8] corresponds to: co=5/3, ci = — 1/12, C2=C3=0. 
The dressed coefficients Cq, Ci are shown in Fig. 3 for A^ = 3. 

4.0 
3.5 
3.0 
2.5 
2.0 
1.5 
1.0 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Fig. 4: Improved coefficients cq and ci (Iwasaki action) 

• The Iwasaki set of parameter values [9] is: co=3.648, ci= — 0.331, C2=C3=0; while, in 
principle, cq and ci depend on go, they are typically kept constant in simulations. The 
corresponding dressed values are plotted in Fig. 4 (A^ = 3). 
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2.6 
2.4 
2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.0 

7.8 7.9 8.0 8.1 8.2 8.3 8.4 8.5 8.6 8.7 

(3 ■ Co 

Fig. 5: Coefficients cq , ci , C3 (red/blue/green dots, respectively) and their dressed counterparts 
Co 5 ci , C3 (dots joined by a line), for different values of /3co = 6co/go (TILW actions) 

• Another class of gluon actions based on Symanzik improvement are the tadpole improved 
Liischer-Weisz (TILW) actions [10,11]. In this case, the coefficients cq , ci , C3 are optimized 
for each value of j3 = 2N / separately (c2 = 0). In Fig. 5 we show the values of q and of 
their dressed counterparts q in a typical range for : 8.0 < (3 cq < 8.6 (A^ = 3). 





Co 


Cl 


Co 


Cl 


1.1636 


5.29078 


-0.53635 


3.39826 


-0.22528 


0.6508 


12.2688 


-1.4086 


8.8070 


-0.7313 



TABLE I. Improved coefficients cq and ci in the DBW2 action, 
for (3 = Q/gl = 1.1636 and 0.6508 



• Finally, the DBW2 gluon action [12] corresponds to C2=C3=0, and /9-dependent values for 
Co, Cl. Some standard values for cq and ci (obtained starting from (3 cq = 6.0 and 6.3), as 
well as Cq and Ci are shown in Table I. 

C. Dressing vertices 

• We will begin by considering the 3-gluon vertex, coming from the action, Eq. (1). This 
vertex results from a Taylor expansion of Ui to 3rd order in g^ . Expressing Ui as in Eq. (4), 
we see that only terms of the form Tr(F/^'* F/^'') will appear in this vertex, since Tr(F/^'') 

and Tr((F/^-') ) will vanish. 

By analogy with Eq. (15), the dressed 3-gluon vertex equals: 
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Consistently with the dressing of propagators, each (21 + l)-point vertex in Eq. (30) is a sum 
of 4 parts (one from each type of Wilson loop in the action) , made up of 

T^iiFl'^f'^F^) (31) 

Denoting the bare 3-gluon vertex by: V3 = Cq V^'^^ + Ci V^*^^'* + C2 v}"^^ + C3 v}^\ it is relatively 
straightforward to see from Eq. (30) that the dressed vertex, V^^, is given by: 

vi'-E-^ ^TlrnrTiv^^ (^^o)-^« (32) 



The summations inside parentheses are a mere multiple of those in Eq. (23); consequently, 
the result for V^'^'' turns out very simple: 

^3'^ = Co V^"^ + ci V^'^ + C2 Vi'^ + C3 yi') (33) 

• We turn now to the 3-point fermion-antifermion-gluon vertex [5] . In the cases of Wilson 
and overlap fermions, these vertices remain unaffected, since the fermion actions do not 
contain any closed Wilson loops on which the BCH formula might be applied. The vertex 
from the clover action, on the other hand, is amenable to improvement; we write: 

^ ■■■ U ■■■ (34) 
(fermions are denoted by a dotted hne). Just as in Eq. (32), we find: 

;gfSi^-(-^^-)^) - -( ■(!) ,35) 



• Vertices with more fields would seem a priori more difficult to handle. To illustrate 
the complications that may arise, let us consider the 4-gluon vertex. The BCH expansion 

of Tr(t/j) contributes to this vertex in the form: Tr((F/^'') ), Tr((F/^'') ) and Ti {f}^^ Fj^^^). 
Such terms may in principle dress differently from each other. In addition, the dressed vertex 

produced from Tr((i^^ ') ) will not be a multiple of its bare counterpart; rather, it will be a 
linear combination of two color tensors (which are independent for N > 3): 

Tr{r"r''r^r'' + permutations} and (S^'^'S^ + S'^'S^'^ + 5'"^5'^) (36) 

This issue has been resolved in Ref. [1], and it generalizes directly to the present case. 

Actually, such complications will not appear while dressing 1- and 2-loop diagrams in typical 

(1) ^ 

cases: Terms of the type Tr((F/ ') ) must simply be omitted in order to avoid double 
counting, since their contribution is already included in dressing diagrams with one less 

loop. Thus, one is left only with: Tr((F/^'*) ) and Tr(F/^'* -F/^'*); for both of these terms it is 
straightforward to show, just as in Eqs. (32, 33), that their dressing amounts to replacing Cj 
by Cj. 

• The same considerations as above apply to all higher vertices from both the gluon and 
fermion actions as well. 
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D. The improvement procedure in a nutshell 



The steps involved in the resummation of cactus diagrams can now be described quite 
succinctly: 

• Substitute gluon propagators in Feynman diagrams by their dressed counterparts. The 
latter are obtained by the replacement Cj — > q = % , where % are the solutions of 



• Perform the same replacement, Ci ^ Ci, on the 3- gluon vertex. 

• Account for dressing of the 3-point vertex from the clover action by adjusting the 
clover coefficient csw '■ csw csw ■ (co/cq) 

• In dressing subleading-order diagrams, avoid double counting, i.e., subtract terms 
which were included in dressing leading-order diagrams. These are very easy to iden- 
tify and subtract: Writing a general subleading-order result (aside from an overall 
prefactor) as: a/N^ + b + cNf/N {Nf : number of fermion flavors), the quantity to 
subtract will include all of a/N"^ (because terms with BCH commutators are higher 
order in N), and it will be a multiple of (2 A?"^ — 3); thus subtraction boils down to the 
substitution: 



(see Refs. [1,13-15] for different applications of this). The remaining subleading ver- 
tices dress exactly as the propagators and 3-point vertices above. 



We turn now to two different applications of cactus improvement: The additive mass 
renormalization for clover fcrmions and the 1-loop renormalization of the axial and vector 
currents using the overlap action. Both cases employ Symanzik improved gluons; hence, our 
results are presented for various sets of Symanzik coefficients. 



It is well known that an ultra-local discretization of the fermion action without doubling 
breaks chirality. Consequently, we must demand a zero renormalized fermion mass, in order 
to ensure chiral symmetry while approaching the continuum limit. For this purpose, the 
bare mass is additively renormalized from its zero tree-level value to a critical value dm. 

We compute dm using clover fermions and Symanzik improved gluons, in 1-loop per- 
turbation theory; the result, denoted as dmi_ioop, is then dressed with cactus diagrams to 
arrive at the improved value c^^^iLioop- The reader can refer to other works [13,14] for more 
details on the definition of dm. A two- loop calculation of dm with the same actions can be 



Eqs. (26). 




(37) 



III. APPLICATIONS 



A. Critical mass of clover fermions 
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found in our Rcf. [15]. There are only two 1-loop diagrams contributing to dmi_ioop and the 
result can be written as a polynomial in the clover parameter: 



2 

(imi_ioop = J^^^'^^sw 

i=0 

^^?-ioop = E4?4w (38) 

i=0 

Clearly, one-loop results are independent of Nf, the number of fermion flavors. e^J^^ includes 
one factor of cq/cq for each power of Csw (cf- Eq (35)). An overall factor of Qq {N'^ — 1)/N 
has been absorbed in the coefficients and e^^^, since the improvement procedure requires 
us to choose definite values of go and N. 

Some numerical values of Eqs. (38) corresponding to the plaquette and Iwasaki actions 
are given below {N — 3). First, for the plaquette action, choosing /3 = 6.0 one gets 



dmi_ioop = -0.43428549(1) +0.1159547570(3) cgw +0.0482553833(1) 4^ 
rfmfii„„p = -0.579221119(2) +0.1159547570(3) cgw +0.0361806779(1) 

in agreement with Ref. [14]. For the Iwasaki action at (5 — 1.95 : 



(39) 



c^^i-ioop = -0.6773690760(3) +0.2342165224(9) cgw +0.0806966864(3) 4^ 
d^i-\oop = -0.757856451(1) +0.1671007819(8) cgw +0.0447467282(1) 4^ 

Eqs. (39, 40) can be used to evaluate the critical hopping parameter through: 



(40) 



i^cr = — r— (41) 

2dm + 8r ^ ' 

where r is the Wilson parameter. Estimates of k^t from numerical simulations exist in the 
literature for the plaquette action [16,17] {Nj = 0), [18,19] {Nj = 2), and also the Iwasaki 
action [20] {Nf = 2). Perturbative (unimproved and dressed) and non-perturbative results 
are listed in Table II for specific values of csw- It is clear that cactus dressing leads to results 
for Kcr which are much closer to values obtained from simulations. 



Action 






Csw 


K 

cr,l — loop 


cr,l— loop 


^non— pert 


Plaquette 





6.00 


1.479 


0.1301 


0.1362 


0.1392 


Plaquette 





6.00 


1.769 


0.1275 


0.1337 


0.1353 


Plaquette 


2 


5.29 


1.9192 


0.1262 


0.1353 


0.1373 


Iwasaki 


2 


1.95 


1.53 


0.1292 


0.1388 


0.1421 



TABLE II. 1-Ioop results and non-perturbative values for k, 
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B. One-loop renormalization of fermionic currents 

As a second application of cactus improvement, we investigate the renormalization con- 
stant Zy [Z a) of the flavor non-singlet vector (axial) current in 1-loop perturbation theory. 
Overlap fermions and Symanzik improved gluons are employed. Bare 1-loop results for Zy^A 
have been computed in the literature [21,6,22]; they depend on the overlap parameter p 

(0 < p < 2). 

One can show that, using the overlap action, the renormalization constants Zy and Za 
are equal [21]; in the MS scheme they read 

C 

^y,A(a, p) = 1 - t/o ^^v,A = 1 - c/o ( + ^t) (42) 

folowing the notation of [6,22]. hy^A and 6^ ai'e 1-loop results pertaining to the amputated 
two-point function of the current and to the fermion self-energy, respectively; since Zy = Za, 
we can write by — Ba- Using cactus improvement, Eq. (42) becomes 

^^>,p)-l-5o'^?y,A (43) 

To compute zfyA dress the Symanzik coefficients and the propagators as described in 
the previous section. In Table III the values of Zy^A and Zyj^ arc presented for different 
sets of the Symanzik coefficients, choosing p = 1.0, p = 1.4. Systematic errors are too small 
to affect any of the digits appearing in the table. The dependence of Zy^A and Zy^^ on the 
overlap parameter p is more clearly shown in Fig. 6, where we plot our results for three 
actions: Plaquette, Iwasaki and TILW. Note that improvement is more apparent for the 
case of the plaquette action. Indeed, from Table III one can clearly see that the effect of 
dressing is smaller for improved gluon actions. This, of course, could have been expected, 
since these actions were constructed in a way as to reduce lattice artifacts, in the ffist place. 



Action 




ZyXp=l-0) 


z^Ap=i.o) 


ZyAp=lA) 


^^a(p=1-4) 


Plaquette 


6.00 


1.26427 


1.35247 


1.14707 


1.19615 


Symanzik 


5.00 


1.24502 


1.29231 


1.13574 


1.16207 


Symanzik 


5.07 


1.24164 


1.28735 


1.13386 


1.15932 


Symanzik 


6.00 


1.20418 


1.23484 


1.11311 


1.13019 


TILW 


3.7120 


1.27581 


1.31941 


1.15259 


1.17690 


TILW 


3.6018 


1.28223 


1.32764 


1.15613 


1.18146 


TILW 


3.4772 


1.28946 


1.33677 


1.16012 


1.18651 


TILW 


3.3985 


1.29434 


1.34298 


1.16282 


1.18995 


TILW 


3.3107 


1.29973 


1.34972 


1.16579 


1.19369 


TILW 


3.2139 


1.30569 


1.35705 


1.16908 


1.19774 


Iwasaki 


1.95 


1.39343 


1.44921 


1.21724 


1.24847 


Iwasaki 


2.20 


1.34872 


1.38773 


1.19256 


1.21440 


Iwasaki 


2.60 


1.29507 


1.31940 


1.16293 


1.17656 


DBW2 


0.6508 


1.49631 


1.45362 


1.27543 


1.25057 



TABLE III. Results for Zy,A,Z^'^^ (Eq. (42,43)), using p=1.0,p=1.4 
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A comparison between our improved Zy^A values and some non-perturbative estimates 
[23], shows that improvement moves in the right direction. Cactus dressing had aheady 
been tested using clover fermions [5], and it turns out to be as good as standard tadpole 
improvement [4], but still not very close to non-perturbative results. 



2.0 




1.0 I 1 1 1 1 ^ ^ 1 

0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

P 

Fig. 6: Plots of Zy^A and Zyj^ for the plaquette, Iwasaki and TILW actions. Labels have been 
placed in the same top-to-bottom order as their corresponding curves. 



In closing, we remark that resummation of cactus diagrams is readily applicable to any 
observable in lattice gauge theories. This procedure for improving bare perturbation theory 
is gauge invariant, and can be applied in a systematic fashion to improve (to all orders) 
results obtained at any given order in perturbation theory. 



APPENDIX A: DRESSING QED 

Cactus improvement can be easily carried over to Lattice Quantum Electrodynamics. In 
this case, link variables commute; hence the first order BCH formula is exact and dressing 
includes the full contribution of diagrams with cactus topology. 
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Dressing the propagator now proceeds precisely as in Eq. (23). The only difference is 
that the result of contracting n—2 out of n generators of SU{N): n F{n; N)/{N'^ — 1), must 

now be replaced by: (n — 3)!! (the number of ways to pair n—2 out of n objects); this 

results in: 

E (n-3)!!2/3r^)/^ = l-e-A^oV2 

n=4,6,8,.. 90 V2/ 

^Ci = Ci-ai = Ci e-^i ^0/2 (Al) 

{Pi are the c-dependent integrals defined in Eqs. (16)). 

As before, the 4 coupled equations (Al) for q assume their simplest form in the case of 
the plaquette action (cq = 1, ci = C2 = C3 = 0); in this case, /3o = l/(2co) and we obtain: 

Co = e-9o/ (4co) ^ ^ ^ Cs ^ (A2) 

Dressing vertices is simpler than in the non-Abelian case. For a bare (2m)-point vertex, 
denoted as: = Cq ^^^^ + Ci + C2 + C3 ¥2^ , instead of contracting group genera- 
tors, one must simply count the number of distinct pairings of 21 objects out of (2Z-|-2m): 

(2/ — 1)!!. The dressed vertex becomes: 

= E^2SQe-A^o/2^f y«g^ (A3) 

1=0 i=0 

The above relations allow us to summarize the dressing procedure for QED very briefly: 

• Replace q by q (as given in Eq. (Al)) throughout 

• Omit all diagrams which contain any cactus subdiagram, to avoid double counting 

Acknowledgements: This work is supported in part by the Research Promotion Foun- 
dation of Cyprus (Proposal Nr: ENTAS/0504/11, ENIEX/0505/45). 
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